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ABSTRACT 

We have observed the cosmic microwave background (CMB) in three regions of sky us- 
ing the Very Small Array (VSA) in an extended configuration with antennas of beamwidth 2° 
at 34 GHz. Combined with data from previous VSA observations using a more compact array 
with larger beamwidth, we measure the power spectrum of the primordial CMB anisotropies 
between angular multipoles ^= 160 - 1400. Such measurements at high £ are vital for break- 
ing degeneracies in parameter estimation from the CMB power spectrum and other cosmo- 
logical data. The power spectrum clearly resolves the first three acoustic peaks, shows the 
expected fall off in power at high £ and starts to constrain the position and height of a fourth 
peak. 

Key words: cosmology:observations - cosmic microwave background 



1 INTRODUCTION 

Acoustic peaks in the power spectrum of cosmic microwave back - 
ground (CMB) anisotrop ies were predicte d by Sakharov (1965), 
ISunvaev & Zel'dovich. tl97ol) and lPeebles & Yu (.igTO) . Recently 
several experiments have accurately measured the first of these 
peaks and have detected the second peak fLee et al j | 200lt 
[Netterfield et al. 2002; Halverson et al. 2002; Scott et ay |20oJ. 
Importantly, despite the differing types of potential systematic er- 
rors suffered by these experiments, the y are in very good agreement 
with each other. The CBI experiment jPadin et aLi(2002|) has mea- 
sured the power spectrum out to £ ~ 3500 ( Mason et alJ l200'S 
^earson et al. 2002), and has detected the predicted fall in power 
level due to incoherent addition of temperatur e fluctuations along 
the line of sight and Silk damping 1968h . but does not have 
the ^-resolution to define the peak structure. Resolving the third 
and subsequent CMB peaks is essential to constrain further cosmo- 



logical models. The Very Small Array (VSA) ('Watson et al.'200j 
(hereafter Paper I) is an interferometric array which has measured 
the CMB power sp ectr um between i = 1 50-800 as described in 
iTavlor et alj <2002l) and lScott et'ZI j2002h (hereafter Paper II and 
Paper III respectively). These measurements were made with the 
VSA's 14 antennas arranged in a compact configuration. The VSA 
has since been upgraded with larger apertures to allow observations 
in an extended configuration giving high temperature-sensitivity 
measurements on smaller angular scales. This paper describes the 
first VSA measurements with the extended array and the power 
spectrum from the combined compact and extended arrays. 

2 THE VSA EXTENDED ARRAY 

As described in Paper I, each VSA antenna comprises a conical cor- 
rugated horn feeding a section of a paraboloidal mirror. In its ex- 
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Table 1. VSA field positions and total effective integration time remaining 
after flagging and filtering of the data. 





RA (J2000) 


DEC (J2000) 


Tint (hrs) 


VSAIE 


00 22 37 


30 16 38 


106 


VSAIF 


00 16 52 


30 24 10 


94 


VSAIG 


00 19 22 


29 16 39 


79 


VSA2E 


09 37 57 


30 41 28 


110 


VSA2F 


09 43 46 


30 41 14 


101 


VSA2G 


09 40 53 


31 46 21 


115 


VSA3E 


15 31 43 


43 49 53 


130 


VSA3F 


15 38 38 


43 50 18 


114 


VSA3G 


15 35 13 


42 45 05 


112 



tended configuration, the 143 mm diameter illuminated apertures 
of the compact array are replaced with 322 mm diameter apertures, 
giving a primary beam of 2.0° FWHM at 34 GHz and an improve- 
ment in flux sensitivity of a factor of just over 5. The filling factor 
of this configuration is greater than that of the compact array, giving 
a significant increase in the overall temperature sensitivity. 

The extended array has several features which increase ob- 
servational efficiency. The narrower primary beam response allows 
observations closer to the Sun and the Moon than for the compact 
array. Repeating the tests outlined in Paper I we conclude observa- 
tions must be at least 18° from the Sun and the Moon. Also the 
higher gain and directivity of the new antennas results in a reduced 
level of cross-coupling, now measured to be less than —120 dB. 
The 'spurious signal' discussed in Paper I has now largely vanished, 
even on the shortest baselines, but conservatively we still apply a 
high-pass fringe-rate filter to our data to remove any low-level con- 
tamination. 

Early commissioning runs with the extended array showed 
that, at the extremes of the VSA's pointing range, the data were 
contaminated when sidelobes of the telescope's primary beam were 
directed towards the top edge of the VSA ground screen; this edge 
appears as a temperature discontinuity and so results in a corre- 
lated signal. Tests show that the contamination is negligible for 
hour angles of less than ±3 hrs; adopting a conservative approach 
we limit observations with the VSA extended array to hour angles 
of ±2.5 hrs. 



3 OBSERVATIONS 
3.1 Field centres 

The observations presented in this paper were made at 34 GHz 
during the period 2001 October - 2002 April. In each of three 
regions of sky, three overlapping fields were observed, separated 
by 75 arcminutes. All the fields lie within the FWHM of the mo- 
saiced fields already observed at lower resolution by the VSA in 
its compact configuration (Tavloretal. 2002). The positions and 
effective integration times of the nine fields are given in Table □ 
The fields were selected for low synchrotron and free-free emis- 
sion, low dust and a lack of large-scale structure and are known 
from previous observations to contain no radio sources brighter 
than 100 mJy at 34 GHz. In Paper II we used external foreground 
templates for synchrotron, free-free and the dust, to estimate the 
level of contamination in the VSA fields. We showed that the am- 
plitude of the combined foreground signal at 1° angular scales is 
< 6 ;UK. It is known that the power spectra of all the Galactic fore- 



grounds decreases with increasing (. with power-law index a = 2-3 
(Q o= r"-) ( Giardino et al. 2001 ; Teamark et al. 2000). Assuming 
a power-law index of 2.5, at £ — 1200 the Galactic foregrounds 
will contribute no more than ~ 16 /jK^ to the power spectrum. We 
therefore conclude that contamination from Galactic emission will 
be negligible compared to our random errors. Foreground contam- 
ination is discussed further by Dickinson et al. (in prep). 

3.2 Data reduction 

The data were reduced and analysed as described in detail in Papers 
I and II but with three differences. 

Calibration. Our primary calibrator was J upiter, whose effec - 
tive temperature at 34 GHz is 734 = 155 ± 5 K ( Mason et al.'l99S(). 
We also use Gas A and Tau A as phase calibrators, transferring the 
flux scale from Jupiter. These sources are partially resolved on the 
longer baselines of the VSA extended array and it is necessary to 
model them. We used V LA 1.4 GHz images (Anderson et al. 199^; 
lBietenholze'tailll99lh with aperture plane coverage encompass- 
ing that of the VSA to model the change in observed flux density 
over the range of VSA baselines. Both Gas A and Tau A are suf- 
ficiently similar in structure at 1 .4 and 34 GHz for this process to 
work well; we have confirmed this using a 32-GHz image from Ef- 
felsberg for Gas A (W. Reich priv comm) and an 850-^(m SGUBA 
map for Tau A jGreerJl2002h . We find that the amplitude of the 
correction is about 15 percent on the longest VSA baselines. 

Effects of bad weather. Some of the extended array data were 
taken in weather conditions considerably worse than those for the 
compact array. Days when the sky temperature increased signifi- 
cantly were identified from the VSA system temperature monitor 
and the data were discarded. However, isolated clouds unresolved 
by the telescope beam can give correlated signals in the visibilities 
without any detectable increase in the system temperature. This 
emission is visible in the data as periods of very rapidly varying 
fringes. These fluctuations are generally short-lived and strongest 
on projected baselines which are perpendicular to the wind direc- 
tion. To detect this correlated emission we calculate the variance of 
the 1 -second sampled data over all 91 baselines during a 16-second 
integration; this is then used to reweight the data for each sample. 
This is in addition to weighting for system temperature variations 
and antenna effective sensitivities. 

Correction for fringe smearing. To remove DC correlator off- 
sets we apply a phase-switching sequence in the hardware with a 
cycle time of 16 s, which is then demodulated in software. For the 
higher fringe rates of the extended array, this process results in a 
reduction in the signal amplitude by a factor of up to 6% and we 
apply an appropriate correction and downweighting of the data for 
each visibility. 

The data were analysed independently by each of the three 
collaborating institutions, and the results were found to be consis- 
tent to within a small fraction of the intrinsic uncertainties in the 
data. In contrast to the compact array analysis, only 28 percent of 
the data were discarded due to weather, fringe rate filtering and tele- 
scope downtime; this difference is due in part to the higher fringe 
rates of the extended array. 

3.3 Source subtraction 

The point source contribution to the observed power spectrum in- 
creases with £ as TqCq£{£+ 1)/2ji, and contamination from extra- 
galactic radio sources is expected to be significant for observations 
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made using the VSA extended array. Our source removal strategy 
is to survey the fields at 15 GHz with the Ryle Telescope (RT) 
and then follow up the sources detected using a single-baseline 
source-subtraction interferometer operating at 34 GHz. The source- 
subtractor baseline is just longer than lOOOA. and so the contamina- 
tion from primordial fluctuations is negligible. 

To ensure that extended-array observations of the CMB are not 
compromised by point-source contamination, it is necessary to sub- 
tract fainter sources than was required for the compact array. Tak- 
ing the 34 GHz source count derived in Paper II, if we subtracted 
sources to the same level as before (i.e. S34 > 80 mjy), then the 
residual source power at £ = 1500 would be Ar5mjj.(,g5 ~ I800/iK^. 
This is significantly higher than the predicted CMB power spec- 
trum. For the current observations we subtract sources down to 
S34 > 20 mJy, which limits the contribution of unsubtracted sources 
to less than ^ AlQfjY? at £ = 1500, and provides better statistics 
for the removal of this residual contribution to the power spectrum. 

To achieve this levl of source removal, several changes were 
made to the source subtraction strategy. The flux limit of the RT sur- 
vey at 15 GHz was improved to 10 mJy, which ensures that sources 
with a spectral index as steep as —1 between 15 and 30 GHz can 
be found (Waldram et al. 2002). The sensitivity of the source sub- 
tractor was also improved by replacing the dishes with ones of in- 
creased surface accuracy and hence higher aperture efficiency, and 
by upgrading to a double sideband system, doubling the effective 
bandwidth. 

In total 24 sources were detected above 20 mJy at 34 GHz 
and removed. We found reasonable upper and lower lim- 
its to the source count in the range 20 - 100 mJy to be 
N{S)dS = 7.6 X 103(5/100mJy)-2 i5jy-lsr-i and N{S)dS = 
8.5 X I0^(5/100mJy)-i '*Jy-isr-' (Taylor et al. in prep. wiU 
present a full discussion of the 34 GHz counts). Integrating these 
counts from zero to 20 mJy gives a range of residual source powers 
of T^Cq = 1.3 ±0.5 X lO^^fjK^. We therefore subtract this statis- 
tical correction from the power spectrum band powers, in addition 
to the subtraction of discrete sources from the visibilities. 



Table 2. Rms noise levels, in mJy beam^', for each of the VSA fields: 
the thermal noise level, Oauto (measured from the autosubtracted maps - see 
Paper II); the noise level in the centre of each map, Cc; the rms CMB signal, 
OcMB- These values are taken from CMB maps made at full resolution. 
The flux-to-temperature conversion is about ImJybeam^' =2.l/jK, giving 
mean rms CMB fluctuations in the maps of Rs 30/jK. 



Field 


^auto 




OcMB 


VSAIE 


11 


19 


15 


VSAIF 


11 


19 


15 


VSAIG 


11 


15 


10 


VSA2E 


10 


19 


15 


VSA2F 


11 


22 


18 


VSA2G 


10 


23 


19 


VSA3E 


10 


18 


14 


VSA3F 


10 


20 


15 


VSA3G 


9 


14 


10 



Table 3. The values for data splits on each of the VSA fields. In each 
case the visibility data from each field were split in two according to epoch 
and the of the difference vectors formed. Also tabulated are the number 
of degrees of freedom (DOF) and the significance of each value. 



Field 


DOF 


X- 


Significance 


VSAIE 


4778 


4851 


0.23 


VSAIF 


3871 


3964 


0.14 


VSAIG 


4070 


3752 


0.99 


VSA2E 


5002 


5170 


0.05 


VSA2F 


3831 


3763 


0.78 


VSA2G 


4314 


4297 


0.58 


VSA3E 


4937 


5011 


0.22 


VSA3F 


4300 


4325 


0.39 


VSA3G 


4970 


5060 


0.18 


Day/Night 


1051 


986 


0.92 



4 RESULTS 

In Figure Q we present maps of the three regions observed. These 
maps provide important data checks and allow comparison with the 
earlier maps of the same sky areas made using the compact array 
but are not used in the calculation of the power spectrum, which 
is computed directly from the visibility data. Individual maps were 

made of each of the nine pointings, using a radial weighting func- 

1/2 

tion in the aperture plane proportional to Cg of a fitted ACDM 
power spectrum - this is equivalent to a Wiener filter. The maps 
were CLEANed to reduce the long-distance correlations arising 
from the restricted sampling in the aperture plane and then com- 
bined into mosaics, weighted by their respective primary beams, 
using the AIPS task Ltess. Figure|2|shows a comparison between 
the compact and extended array maps of the VSAl field; the two 
maps have different but overlapping aperture plane coverages and 
show many features in common. 

We find that the noise levels in our maps well outside the pri- 
mary beam are consistent with the thermal noise level derived from 
subtracting alternate visibilities (CTauto)^ showing that there is no 
component in the maps other than the signal within the primary 
beam and the thermal noise. The thermal noises for each field are 
given in Table |2| together with the rms CMB contributions in the 
central region of each of the maps. 



4.1 Data checks 

As a check on the consistency of the data, we computed the 
statistic and significance for splits on the visibility data following 
the method described in Paper III. The significance is given as the 
probability to exceed the observed value in the cumulative dis- 
tribution function. The data for each VSA field were split in two 
according to observing epoch and the values and associated sig- 
nificances are given in Table|3| A split between day and night ob- 
servations is also included. 

The consistency of the power spectra derived from each of the 
3 VSA mosaiced regions was compared by forming the statistic 
on pairs of power spectra. The x^ values (and significances) for the 
VSA1/VSA2, VSA1/VSA3 and VSA2A^SA3 power spectra com- 
parisons are 12.2 (0.73), 9.5 (0.89) and 15.4 (0.50) respectively In 
each case there are 16 degrees of freedom in the power spectrum 
analysis. The data sets are clearly self-consistent. 



5 POWER SPECTRUM 

The fully reduced and source subtracted visibility data for the com- 
bined compact and extended VSA data sets have been analysed us- 
ing the Madcow software package I Hobson & Maisinaer 200^. 
The typical anti-correlation between adjacent bins is ~ 20 percent. 
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Figure 1. Mosaiced and source-subtracted maps for each of the VSA regions: left, VSAl; centre, VSA2; right, VSA3. In each case 
the greyscale runs from —150 to 150/jK and the contour indicates the half-power response of the combined primary beam. 




Figure 2. Comparison of the compact and extended array maps for the cen- 
tral part of the VSAl field. Greyscale: Compact array map, with greyscale 
running from —250 to 250^K. The sensitivity is nearly uniform across this 
area. Contours: Extended array map, contour interval 30/.iK. The map has 
been corrected for the combined effect of the primary beams, so the noise 
level rises towards the edge of the map. 



with one as high as 40 percent. We use variable sized bins to reflect 
the differences in density of coverage of the aperture plane, but also 
repeat the calculation with bin centres shifted by half a bin width in 
order to sample effectively the features of the power spectrum. We 
follow the procedure outlined in Paper III to determine the window 
function W{t} which determines how each bin samples the under- 
lying power spectrum. The window functions for our analysis are 
shown in Figure|3| The resulting power spectrum is shown in Fig- 
ure|4| numerical values for both binnings are given in Tableland 
the correlation matrix for the main binning in Table|5] The bin cen- 
tres given in each case are the medians of the respective window 
functions rather than the nominal bin centres. 

There is a systematic uncertainty of 7 percent in the scaling 
of the power spectrum due to the uncertainty in the effective tem- 
perature of the primary calibration source, Jupiter. The error in flux 
density to temperature conversion (equivalent to the beam uncer- 
tainty in scanned-beam experiments) is negligible. 




Figure 3. Window functions for the combined compact and extended array 
data set. The functions are normalised to unit area, and different bins are 
plotted with different Unestyles. 



6 DISCUSSION & CONCLUSIONS 

We have measured the power spectrum of the CMB over the mul- 
tipole range I = 160 - 1400 using two configurations of the VSA. 
The combined power spectrum clearly shows three peaks, a sharp 
fall-off in power above £ ~ 800, and marginal evidence for a 
fourth peak. Figure |5| shows the VSA data plotted alongside that 
from other recent CMB experiments: there is excellent agreement 
between all the experiments. 

To quantify the significance of these peaks in the power spec- 
trum, we used a model of five Gaussians parametrised by their am- 
plitude, position and width, as described in Odman et al. (200^. 
We analysed the VSA data alone, plus two data sets consisting of a 
compilation of recent data, and used an MCMC routine to constrain 
the parameters simultaneously. In Table|5| we list the constraints on 
the amplitudes and positions of the first four peaks after marginal- 
ising over their widths (there are no interesting constraints on the 
fifth peak). Including the new VSA measurement improves the con- 
straints on the amplitude of the fourth peak, as well as confirming 
the detection of the first three peaks. 

Fits of adiabatic inflationary model power spectrum, with as- 
sociated estimates of cosmological p arameters, are presented in a 
companion paper iSlosar et all2002fl . 

In our current observing configuration we are increasing the 
size of the mosaiced regions and also observing new regions; this 
will increase both the sensitivity and i'-resolution of the power spec- 
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B 



^'-ranee 



4 



1 

lA 
2 

2A 
3 

3A 
4 

4A 
5 

5A 
6 

6A 
7 

7A 
8 

8A 
9 

9A 
10 

lOA 
11 

llA 

12 
12A 

13 
13A 

14 
14A 

15 
15A 

16 



100- 
145- 
190- 
220- 
250- 
280- 
310- 
340- 
370- 
410- 
450- 
475- 
500- 
540- 
580- 
610- 
640- 
670- 
700- 
725- 
750- 
800- 
850- 
900- 
950- 
1000- 
1050- 
1125- 
1200- 
1275- 
1350- 



190 
220 
250 
280 
310 
340 
370 
410 
450 
475 
500 
540 
580 
610 
640 
670 
700 

-725 
750 
800 
850 
900 
950 

1000 

1050 
1125 
1200 

- 1275 
1350 

- 1525 
1700 



160 
190 
220 
251 
289 
321 
349 
376 
416 
431 
479 
501 
537 
581 
605 
639 
670 
696 
726 
759 
795 
843 
888 
948 
1002 
1057 
1119 
1199 
1271 
1357 
1419 



T^Ht + l)Ctl2n\ftY.^] 



3864115^1 
5893llf37 
5390+'-'*'' 
2603li« 
17491^47 

16381^^^ 
286615« 
1460lg^ 
22371^^ 



19221™! 

35871^^1 

147llg^ 
n+1091 

1 1 25+*''* 



(,+ 1431 



1311+'"* 



7031 + '^*' 
3450l?« 
19541^^^ 
15971^1^ 
2675lg^ 
22101^^ 
15011'^^ 
23281™ 
29051^^1 
28761^^ 

20715;;^ 
509IJ4I 

5901^^3 
505li^r 



Table 6. Constraints on parameters of a multiple Gaussian fit to the 
power spectrum data. The inclusion of the new VSA data places con- 
straints on the power beyond the first three peaks. Dataset I consists of 
the QMASK compilation (Xu et al. 2002), MSAM fwilson et al. 20q3), 
Python V 1 Coble et al. 2001), Boomerang 98 ( Netterfield et al. 2002), Max- 



ima iL eeetal,, 2001.). DASI i Haly erson et al. 2002) , CBI 1 Pearso n et alj 
2002), ArcheoDs TBenoit et all2002l) and the VSA data from the compact 
aiTay configuration. Dataset II is identical, except for the inclusion of the 
new VSA data. Note that, while the lack of a lower bound on a AT implies a 
significant contribution to the power from adjacent Gaussians, naiTow con- 
straints on peak locations imply strong evidence for spectral features rather 
than a flat spectrum of excess power. 



VSA only 



dataset I 



dataset II 







70.2li:i 


71.7+I'? 




45.5li^ 


38.5+^-^ 


42.3lto'4 


AT3 


47 0+** '' 


42.5+[^ 


43.31^] 




25 9+"'-' 

^unbounded 


IS n+" ' 


19.6+||'^ 




224I2J 


214+^ 


2161^ 




5091*1 


51111^ 


51011' 


h 


771+55 


7861^^ 


78211 




1171 +*' 


12601^1 


115711?^ 



trum. We also plan to increase the ^-range of the VSA using larger 
antennas and longer baselines to study further peaks in the primor- 
dial spectrum and giv e information o n the region where CBI ob- 
serves excess power (Mason et al. 2002). 

Band powers, correlation matrices and window functions can 
be downloaded from 

http : //mrao . cam. ac .uk/teles copes/ vsa/results .html 



Table 4. The bandpowers for main and offset binnings of the complete VSA 
data set, combining both compact and extended array data. The <!-range 
gives the nominal bin hmits, while is the median value of the relevant 
window function. The error bars enclose 68 percent of the likelihood, even 
when the lower limit is zero. 



B 


Cb.B-2 


Cb.b- 1 


Cb.b 


Cb.b+ 1 


Cb.B+2 


Co\B.B 


1 






1.0 


-0.09 


0.012 


17.039 


2 




-0.09 


1.0 


-0.206 


0.021 


22.295 


3 


0.012 


-0.206 


1.0 


-0.114 


0.021 


14.886 


4 


0.021 


-0.114 


1.0 


-0.151 


0.031 


3.442 


5 


0.021 


-0.151 


1.0 


-0.227 


0.044 


1.316 


6 


0.031 


-0.227 


1.0 


-0.254 


0.047 


3.584 


7 


0.044 


-0.254 


1.0 


-0.268 


0.037 


3.912 


8 


0.047 


-0.268 


1.0 


-0.211 


0.09 


4.011 


9 


0.037 


-0.211 


1.0 


-0.406 


0.092 


5.135 


10 


0.09 


-0.406 


1.0 


-0.286 


0.055 


7.863 


11 


0.092 


-0.286 


1.0 


-0.187 


0.049 


4.883 


12 


0.055 


-0.187 


1.0 


-0.251 


0.079 


3.905 


13 


0.049 


-0.251 


1.0 


-0.336 


0.052 


2.353 


14 


0.079 


-0.336 


1.0 


-0.162 


0.03 


5.475 


15 


0.052 


-0.162 


1.0 


-0.196 




6.714 


16 


0.03 


-0.196 


1.0 






27.092 



Table 5. The correlation matrix Ci. j for the combined VSA data (main bin- 
ning only). The values of the matrix for which Cj.j is not reported can be 
assumed to be zero. The final column gives the diagonal elements of the 
covariance matrix in units of 10^ x fjK* 
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Figure 4. Combined CMB power spectrum from the three mosaiced VSA fields. Tiie error-bars represent la limits; the two sets of data points correspond to 
alternative interleaved binnings of the data. 




Figure 5. A comparison of the VSA data (filled circles) with results from the BOOMERANG (open squares), MAXIMA (open circles), DASl (open triangles) 
and CBI (stars) experiments. Two sets of error bars are plotted for each data set; the smaller (inner) of the two indicate only random errors, whilst the larger 
bars indicate the amount by which the itmer points could move due to absolute calibration and beam uncertainty. In each case the error bars indicate la limits. 
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